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ABSTRACT 

The  experimental  apparatus  for  temperature  programmed  thermogravimetry 
has  been  modified  to  more  effectively  obtain  kinetic  parameters  for  the 
degradation  of  polymers.  The  thermobalance  was  modified  to  incorporate 
direct  sample  temperature  measurement  thereby  to  minimize  temperature 
measurement  errors.  An  automatic  data  acquisition  system  was  incorpo¬ 
rated  into  the  apparatus  and  appropriate  computer  programs  to  handle  the 
magnetic  tape  data  were  written.  The  modified  apparatus  has  been  tested 
with  several  polymer  systems  and  it  was  demonstrated  that  the  use  of  the 
magnetic  tape  data  recording  system  permitted  greatly  increased  output 
from  the  thermobalance. 
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SECTION  I 
INTRODUCTION 


In  the  previous  report  (Reference  1)  a  method  of  obtaining  kinetic 
parameters  for  the  degradation  of  polymers  using  temperature  programmed 
thermogravimetry  was  described.  The  experimental  procedures  and  a  method 
of  processing  TGA  data  on  the  computer  were  described  including  the 
application  of  the  technique  to  several  polymers .  The  technique  has 
since  been  applied  to  a  variety  of  polymer  systems  with  considerable 
success  (References  2  and  3) .  Routine  operation  of  the  system  revealed 
two  possible  limitations  to  the  accuracy  and  usefulness  of  the  apparatus 
in  its  present  form. 

1.  The  temperature  of  the  degrading  sample  was  assumed  to  be  that 
of  a  thermocouple  placed  near  the  crucible  with  some  temperature 
correction  applied. 

2.  The  output  was  limited  by  the  speed  at  which  data  could  be  read 
off  the  chart  and  prepared  for  processing  by  the  computer. 

Since  the  system  had  been  shown  to  be  capable  of  producing  high 
quality  data,  it  seemed  desirable  to  redesign  the  apparatus  to  remove 
these  limitations  on  its  use.  This  is  described  in  detail  in  the 
following  sections. 
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SECTION  II 

MODIFICATION  TO  THE  AINSWORTH  RV  THERMOBALANCE  TO 
INCORPORATE  DIRECT  SAMPLE  TEMPERATURE  MEASUREMENT 

1.  INTRODUCTION 

In  TGA  it  is  customary  to  calibrate  the  temperature  inside  the  sample 
holder  against  an  external  thermocouple  placed  as  close  to  the  operating 
position  of  the  sample  holder  as  possible,  under  normal  run  conditions 
(heating  rate,  etc.)  except  that  weight  is  not  being  recorded.  Providing 
the  same  conditions  are  observed  during  the  normal  run  there  is  no 
reason  to  suppose  this  technique  is  inaccurate.  However,  for  a  large 
number  of  samples,  heating  rates,  etc.,  this  represents  an  inordinately 
large  number  of  calibrations  and  this  still  presupposes  absolute 
reproducibility  of  the  two  runs.  A  much  more  satisfactory  method  is  to 
measure  the  temperature  of  the  sample  directly  during  the  degradation, 
particularly  in  kinetic  studies  where  temperature  is  so  important.  The 
evaporation  of  material  from  degrading  polymers  can  cause  considerable 
decrease  in  sample  temperature,  particularly  when  rate  of  weight  loss  is 
high.  For  example,  polytetraf luoroethylene  loses  16%/minute  at  its 
maximum  rate  of  weight  loss  under  the  conditions  used  to  study  this 
polymer . 

2 .  MODIFICATIONS 

The  Ainsworth  RV  thermobalance  used  in  this  work  is  particularly 
suitable  for  conversion  to  direct  sample  temperature  measurement.  All 
of  the  parts  are  accessible  when  the  cover  is  removed.  The  fact  that 
the  balance  is  not  the  null  deflection  type  poses  some  problem  since,  at 
some  stage,  wires  have  to  be  taken  from  the  beam  to  a  measuring  device 
thereby  interfering  with  the  normal  free  swing  of  the  balance.  The 
configuration  of  the  wires  described  in  this  section  was  arrived  at  by 
trial  and  error. 

Figure  1  shows  a  general  view  of  the  balance  with  the  bell  jar  in 
place.  Figure  2  shows  the  detailed  arrangement  of  the  connecting  wires. 
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Figure  2.  Details  of  Wires  Attached  to  Balance  Beam  and  the 
Thermocouple/Support  Arrangement 
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Originally  the  balance  was  wired  with  Chromel/Alumel  wire  because  of  its 
high  millivolt  per  degree  output  (0.04  mv/°C) ,  but  because  the  wire  is 
magnetic  there  was  considerable  interaction  with  the  furnace  electrical 
supply  up  to  300°C.  Weight  readings  below  this  temperature  could  not  be 
used.  The  balance  was  later  rewired  with  platinum-platinum/10%  rhodium 
wires  which  are  nonmagnetic.  The  lower  EMF  was  measured  with  a  Digital 
Voltmeter.  The  real  problem  in  direct  sample  temperature  and  weight 
measurement  is  in  finding  some  way  of  transferring  the  EMF  signal  from 
the  balance  beam  without  interfering  with  the  weighing  characteristics. 
Any  attachments  to  the  beam  have  the  potential  of  upsetting  both  the 
sensitivity  and  the  zero  of  the  balance,  especially  as  the  Ainsworth  is 
not  a  null  deflection  balance.  The  wire  attachments  are  shown  in 
Figure  2.  Loop  1  joins  the  thermocouple/suspension  to  the  beam. 

Consider  the  effect  of  changes  in  sample  weight  on  this  loop.  The 
situation  is  shown  in  Figure  4.  <2 1 ,  is  the  angle  between  the  beam  and 

suspension,  initially,  and  a ^  the  angle  after  the  sample  has  lost 
weight.  It  can  be  seen  that  the  arrangement  of  Loop  1  will  tend  to 
restrict  this  motion  and  cause  anomalous  weight  readings.  The 
arrangement  of  Loops  2  and  3  will  have  a  similar  effect.  Loop  3  was 
found  to  have  a  profound  effect  on  the  zero  of  the  balance.  Careful 
arrangement  of  the  length  and  position  of  the  wires  resulted  in  a  stable 
system  provided  certain  limitations  were  recognized.  Although  this 
balance  can  follow  weight  changes  up  to  200  mg  using  multiple  chart  scans 
the  weight  loss  that  could  be  followed  was  less  than  10  mg,  i.e.  one  span 
of  the  recorder  chart.  Some  relaxation  effects  were  noted  when  the  beam 
was  switched  from  one  position  to  another  by  adding  or  subtracting  10  mg 
to  the  counter  weight.  This  limitation  of  10  mg  samples  is  of  little 
importance  since  heating  effects  and  diffusion  usually  restrict  sample 
size.  A  detailed  drawing  of  the  support  suspension  is  given  in  Figure  3. 
Removing  and  rehanging  the  suspension  was  found  to  give  small  variations 
in  readings  so  the  balance  could  not  be  used  to  measure  absolute  sample 
weight.  In  practice  this  was  not  a  problem  since  either  the  sample 
degraded  completely  or  the  weight  added  to  the  crucible  could  be 
measured  with  sufficient  accuracy. 
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The  operation  of  the  balance  is  very  convenient  for  checking  zero 
and  sensitivity,  and  long-term  drift  in  both.  If  the  balance  is  adjusted 
to  give  a  zero  or  100%  reading  on  the  chart,  adding  or  subtracting  10  mg 
by  the  remote  control  switches  the  balance  from  one  extreme  to  the  other. 
In  this  way  changes  in  zero  and  sensitivity  can  be  detected  and  adjusted. 
In  practice,  as  well  as  checking  the  sensitivity  before  a  run,  the 
weight  was  arranged  such  that,  at  the  end  of  the  degradation,  the  balance 
was  close  enough  to  the  zero  position  to  allow  switching  and  a  further 
check.  Sensitivity  variations  were  usually  less  than  1%.  The  constancy 
of  buoyancy  correction  is  another  "built-in"  check  on  the  accuracy  of 
operation. 

The  weight  of  the  Teflon  block  and  quartz  rod  in  the  suspension  was 
found  to  give  good  electrical  contact  at  the  hooks.  The  black  wax  used 
as  an  anchor  for  the  wires  on  the  balance  beam  also  acted  as  an 
effective  insulator.  To  check  the  electrical  integrity  of  the  system, 
the  sample  temperature  as  measured  by  the  suspension  was  checked  against 
an  independent  thermocouple  in  the  crucible.  Variation  was  less  than 
1°C  at  600°C . 

3.  TESTING  UNDER  RUN  CONDITIONS 

Since  an  extensive  study  of  the  degradation  of  polytetraf luoro- 
ethylene  had  been  made  on  the  unmodified  balance  (Reference  1)  a 
complete  kinetic  analysis  was  carried  out  on  the  polymer  using  the 
balance  with  the  wires  attached. 

Samples  of  Teflon  molding  powder  (8-9  mgs)  were  degraded  at  nominal 
heating  rates  of  75,  150,  300,  and  450°C/hour.  The  data  was  analyzed 
by  the  standard  procedures  detailed  in  Reference  1. 

A  plot  of  Activation  Energy  against  %  weight  loss  is  shown  in 
Figure  5,  before  and  after  balance  modifications.  The  results  of  this 
work  show  an  average  activation  energy  of  59.5  kcal  compared  with  69.3 
kcal  for  the  previous  work  (both  for  10-80%  of  the  reaction) .  The 
earlier  results,  however,  show  a  considerable  increase  in  activation 
energy  after  50%  weight  loss. 
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Figure  6.  Arrhenius  Plots  for  PTFE  at  50%  Weight  Loss 
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A  better  comparison  of  the  result  is  obtained  by  comparing  the 
Arrhenius  plots  for  the  two  series  at  50%  weight  loss.  The  points  at 
lower  heating  rates  in  both  cases  fall  on  parallel  straight  lines 
indicating  the  same  activation  energy.  The  separation  of  the  lines 
represent  a  temperature  difference  of  13°C,  the  sort  of  difference  one 
might  expect  between  thermocouples  placed  in  and  adjacent  to  the  sample. 
It  is  interesting  to  note  that,  with  the  thermocouple  in  the  sample,  the 
data  at  300°C/hour  heating  rate  falls  on  the  straight  line  whereas  it 
does  not  in  the  previous  data.  This  is  probably  due  to  the  temperature 
in  the  sample  being  lower  than  that  recorded  in  the  earlier  work. 
Deviations  occur  in  both  cases  at  450°C/hour.  At  this  heating  rate  the 
rate  of  volatilization  is  of  the  order  of  16%  per  minute  and  questions 
of  how  well  the  thermocouple  can  respond  to  the  changes  and  how  the 
sample  is  distributed  with  respect  to  the  thermocouple  arise.  There 
also  exists  the  possibility  of  lower  rates  due  to  diffusion  effects  at 
high  heating  rates,  particularly  in  the  larger  samples  used  in  the 
earlier  work  (100  mg) .  This  effect  of  sample  size  may  have  something  to 
do  with  the  otherwise  unexplained  increase  in  activation  energy  after 
50%  reaction,  observed  in  the  earlier  work. 

4 .  CONCLUSION 

In  general  the  agreement  in  the  two  sets  of  data  is  good  indicating 
that  the  attachments  to  the  balance  beam  have  had  little  effect  on  the 
accuracy  of  the  system.  The  modified  system  is,  however,  inherently 
more  accurate  since  the  temperature  sensor  is  inside  the  crucible, 
although  question  may  still  arise  about  contact  with  the  sample,  thermal 
conductivity  of  the  sample,  temperature  gradients  and  heat  being 
conducted  away  from  the  sample  by  the  wires  (Reference  4) . 

Confidence  in  the  stability  and  response  of  the  balance  was  further 
increased  when  a  set  of  data  obtained  with  a  chromel-alumel  system  with 
manual  reading  of  data  from  a  recorder  chart,  gave  the  same  kinetic 
parameters  for  BBB  degradation  as  the  same  balance  wired  with  platinum- 
platinum/10%  rhodium  and  using  a  magnetic  tape  recording  of  the  data. 
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SECTION  III 

COLLECTION,  PROCESSING  AND  ANALYSIS  OF  TGA  DATA 

1.  INTRODUCTION 

In  the  previous  reports  (References  1  through  3)  TGA  data  was 
obtained  by  reading  several  hundred  sets  of  weight/temperature  data 
points  from  the  recorder  chart,  and  having  the  data  transferred  to 
punched  cards  for  processing  by  the  computer.  This  operation  was  both 
time  consuming  and  tedious  and  considerably  reduced  the  amount  of  data 
that  could  be  produced.  The  method  had  also  considerable  potential  for 
human  error.  Modern  advances  in  instrumentation  have  suggested  the 
replacement  of  the  chart  recorder  by  another  device  such  as  a  magnetic 
tape  or  paper  tape  recorder  which  can  be  read  directly  by  the  computer. 
For  this  purpose  an  SRL  Model  837  Data  Aquisition  System  was  acquired. 
This  is  described  in  the  next  section. 

2.  THE  SRL  MODEL  837  DATA  AQUISITION  SYSTEM 

A  block  diagram  of  the  apparatus  is  shown  in  Figure  7.  The  system 
consists  of  the  following  components: 

1.  Two  model  2670  Data  Amplifiers  -  Hewlett-Packard. 

2.  A  model  X-2P  Digital  Voltmeter  -  Non-Linear  Systems,  Inc. 

3.  A  model  1600  Incremental  Tape  Recorder  -  Kennedy. 

4.  Scanner  and  Counter  Logic  -  SRL  design  using  Digital  Equipment 
Corporation  Flip  Chip  Modules  and  power  supply. 

The  complete  system  is  housed  in  a  67-inch  Honeywell  modu-mount 
enclosure  and  each  basic  component  has  its  own  power  supply,  switch, 
and  fuse. 

Electrical  signals  proportional  to  the  weight  and  temperature  are 
fed  to  the  two  Data  Amplifiers,  the  levels  of  which  can  be  adjusted  to 
send  a  measurable  output  to  the  digital  voltmeter.  The  two  signals  are 
scanned  alternately,  the  scanning  interval  varying  from  0.5  to  10 
seconds  (i.e.  the  interval  between  two  successive  weight  or  temperature 
readings  can  be  varied  from  1  to  20  seconds) .  If  necessary  a  permanent 
record  of  the  data  can  be  obtained  from  the  printer  in  which  case  the 
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lower  limit  of  the  scan  interval  is  governed  by  the  tracking  speed  of  the 
printer.  In  practice  the  printer  is  only  used  during  testing  or 
trouble-shooting.  The  amplified  signals  are  fed  to  the  magnetic  tape 
recorder.  Data  is  recorded  in  records  of  a  length  determined  by  the 
control  logic,  with  a  record  gap  at  the  end  of  each  record.  At  present 
the  apparatus  is  set  up  to  receive  18  sets  of  weight/temperature  data 
but  this  can  be  varied.  This  short  length  is  very  suitable  for 
correction  as  will  be  discussed  later.  When  the  "Stop"  button  of  the 
system  is  activated,  recording  continues  to  the  end  of  the  record.  At 
the  end  of  the  last  record  an  "End  of  File"  code  must  be  recorded.  This 
is  used  by  the  computer  to  detect  the  end  of  the  data  and  without  it 
data  cannot  be  recovered . 

Once  the  data  is  recorded  it  is  now  in  a  form  suitable  for  processing 
on  the  IBM  7094  computer. 

To  minimize  the  loss  of  data  which  could  occur  due  to  various 
failures,  each  run  is  recorded  on  a  separate  magnetic  tape  (Ampex  Data 
Mailer,  200  ft).  Since  it  is  necessary  to  retain  data  for  some  time  but 
undesirable  to  accumulate  numerous  magnetic  tapes,  the  data  is 
transferred  to  a  master  storage  tape  during  the  initial  processing  A 
block  diagram  of  the  tape  manipulation  is  shown  in  Figure  8, 


Figure  8.  Transfer  of  Data  From  Small  Tape  to  the  Master  Storage  Tape 
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Transfer  and  storage  of  data  is  carried  out  when  the  data  is  sent 
for  preliminary  examination  using  Program  1  (Appendix) .  To  minimize 
loss  of  data  due  to  machine  or  operation  error,  three  master  tapes  and 
at  least  the  six  most  recent  runs  are  retained.  The  three  storage  tapes 
have  n,  n-1  and  n-2  runs.  When  hung  in  the  configuration  shown,  the 
tape  with  the  largest  number  of  runs  is  the  old  master  tape  which  is 
read  only.  The  tape  with  the  least  number  is  the  new  master.  The  n-1 
tape  is  meanwhile  safely  stored.  The  data  from  the  old  master  is 
written  on  the  new  master  (Step  1)  followed  by  the  current  run  data 
(Step  2).  This  tape  then  becomes  the  main  master  storage.  The  method 
also  allows  the  erasure  of  the  latest  record  should  the  output  show  the 
data  was  unsatisfactory. 

As  well  as  handling  the  storage  of  data.  Program  1  also  displays  the 
weight  loss,  temperature,  and  rate  of  weight  loss  at  each  of  the  weight 
losses,  along  with  a  record  by  record  account  of  the  data  as  stored  on 
the  tape.  Both  the  records  and  the  number  of  data  points  are  counted 
and  those  figures  are  particularly  useful  for  identification  purposes  in 
the  case  of  bad  data.  Three  types  of  bad  data  have  been  encountered  and 
Program  1  is  available  with  modification  to  cope  with  each: 

1.  Redundant  records  caused  by  write  errors,  eg  parity  errors, 

in  the  recording.  Provided  they  do  not  occur  at  critical  stages  in  the 
degradation  up  to  nine  records  can  be  discarded.  This  type  of  failure 
is  recognized  by  the  computer  in  reading  the  tape  and  the  number  of 
redundant  records  is  shown  on  the  initial  print  out. 

2.  Records  which  have  bad  data  but  which  are  not  redundant  and  are 
not  detected  by  the  machine.  If  they  do  not  occur  at  a  critical  stage 
in  the  degradation  they  can  be  discarded. 

3.  Bad  data  points  in  a  record.  These  can  be  replaced  by  values  in 
keeping  with  the  rest  of  the  data. 

Once  a  set  of  satisfactory  runs  have  been  loaded  on  to  the  master 
tape,  the  data  is  reprocessed  using  Program  2.  This  program,  provides  a 
print  out  of  the  rate  of  weight  loss  at  1%  intervals  and  also  gives  the 
output  on  IBM  punched  cards  for  use  in  the  Arrhenius  Program  (Program 
3 ,  Appendix) . 
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3.  CONCLUSION 

The  complete  series  of  modifications  to  the  thermogravimetric  system 
described  in  this  report  have  been  tested  on  a  series  of  styrene-acrylo¬ 
nitrile  copolymers.  These  copolymers  have  a  very  high  rate  of  weight 
loss  providing  an  effective  test  for  the  direct  sample  temperature 
measurement.  As  with  the  degradation  of  Teflon,  it  was  shown  that  good 

Arrhenius  plots  could  be  obtained  provided  the  heating  rate  did  not 
o 

exceed  300  C/hour.  The  short  degradation  time  was  useful  in  testing  the 
efficiency  of  the  data  recording  system.  It  was  clearly  demonstrated 
that  the  use  of  the  magnetic  tape  data  recording  system  permitted 
maximum  output  from  the  thermobalance.  The  system  is  set  up  such  that 
if  the  output  is  to  be  further  increased  a  second  thermobalance  could  be 
readily  accommodated,  one  balance  being  loaded  and  evacuated  while  the 
other  is  being  used.  Detailed  results  of  the  analyses  of  the  kinetics 
of  degradation  of  the  styrene-acrylonitrile  copolymers  will  be  described 
in  another  report. 


16 


AFML-TR-68-181 
Part  III 

REFERENCES 

1.  I.  J.  Goldfarb,  A.  C.  Meeks,  R.  McGuchan,  AFML-TR-68-181,  Part  II. 

2.  I.  J.  Goldfarb,  R.  McGuchan,  AFML-TR-68-182,  Parts  I  and  II. 

3.  I.  J.  Goldfarb,  A.  C.  Meeks,  AFML-TR-68-367 ,  Part  I. 

4.  H.  C.  Anderson,  in  Thermal  Analysis,  Vol.  1,  Ed.  P.  E.  Slade  and 

L.  T.  Jenkins. 


17 


AFML-TR-68-181 
Part  III 


APPENDIX 


PROGRAMS  1-3 
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PROGRAM  1 


This  program  transfers  current  run  data  to  master 
storage  tape,  provides  record  by  record  output  of 
data  as  it  appears  on  the  tape,  and  provides  a 
preliminary  print  out  of  the  rates  for  examination. 
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TGTAPE  09/17/7C 

_ T  G  T  A  P  E_ JE.EM .  SOU  RCE  ST  ATE  MINT  IF  N  (  S )  - _  _ _ 


C  PROGRAM  1  TAPE  PREPARATION  AND  INITIAL  EXAMINATION  OF  DATA 

C _ PROGRAM.  J.C  RJE AC  CURRENT  DATA.  AND  TRANSFER  IT  JO  THE  MASTER  TAPE 

C  ALONG  WITH  ALL  THE  DATA  ON  THE  OLD  MASTER 

__£  .  THE  PRCORAM..ALSO  PRINTS .  CUT.  THE  CAT  A  OF  THE  CURRENT  RUN  AND 

C  DETERMINES  RATE  OF  WEIGHT  LOSS  AT  ONE  PER  CENT  INTERVALS 

C .  UNIT  ),...=  SM  A  |  L.  TAPE  WI JK.  f.Ufi&ENI,  SUM . . . . .  . . ______ . 

C  UNIT  2  =  NEW  MASTER  TAPE 

...G. .  .  -UNIT  3  ?  OLE  MASTER  ..T APE  (CAN  BE  READ_QNLY  ).  .  ._ _ _  _ 

C  PROGRAM  USES  SUBROUTINE  EOF  TO  PERMIT  READING  OF  A  NUMBER  OF 

C _ FILES  SEPARATE  CL  BY  .  E_N  0  OF  FILE  MARKERS _ _ _ . 

C  INPUT  TEMPERATURES  ARE  FITTED  TO  A  FIFTH  DEGREE  POLYNOMIAL  USING  A 

_C _ LEA-SI — SG.UAR  E  S  SUBROUTINE  (PESO). _ _ _ __________  . _ _ ______ _ 

C  WEIGHTS  CORRESPONDING  TO  SHORT  TEMPERATURE  RANGES  ARE  FITTED  TO  A 

_C _ CUACRATIC  EY  PESO. _ _ _ _  ..  ..  _ 

C  INPUT  WEIGHTS  DIFFERING  FROM  FITTED  LINE  BY  MORE  THAN  ONE  PERCENT  OF  THE 

C TOTAL  WEIGHT  LOSS  ARE  REPLACED  BV  THE  CURVE  FIT  VALUE. _  _  __ 

C  W  =  WEIGHT  CATA  POINT  READ  OFF  TAPE 

r.  _ t  =  tfmpfratiirf  data  point  »i  Ar.  ff  tape i  in  mv.  i _ , _ _________ _ _ 

C  T I  =T IME  DATA  POINT  CALCULATED  FROM  TIME  INTERVAL  AND  NO  OF  CATA  POINTS 

_C  ...  .  T.I  N.T*.tJlSE_  IN  T.ERV A  L _ _ _ _ _ _ _ 

C  WW  =WE IGHT  LOSS  AT  1  PER  CENT  INTERVALS 

C  ...  CWCT=RATE  CF  WEIGHT  LOSS  AT  L  PER.  CENT  INTERVALS  . . .  _  _ 

C  TOE  R=HE  AT  I NG  RATE  AT  1  PER  CENT  INTERVALS 

-C - IPCLY-.LEMPERAT.U.R.E.  IM ESE.QMDIMG._T.O..^AnH- J^ER. _CEN T  WEIGHT  _LOSSf CALCULATED  BY 

C  PLSQ 

...£._ . I N  W  = ._  T.I ME.I1C R.R£S.P.(1N D ING._T.0j.  EACH  PER.  CENT  HEI.6HI. .U.QS.L _ _ 

C  PLOT  =D  IMENS IONS  FCR  GRAPH  PLOT  SUBROUTINE 

C .  B  =  C.G.EFF  10. 1  ENTS  OF  10TH  ORDER  POLYNOMIAL  FITTING  TEMP/EMF  DATA  FOR 

C  PLATINUM*PLATINLM  1CPER  CENT  RHODIUM 

_£. - C ..?.CG.E£f J..CI.E.NJLS.  .OF  WE1GKI/T1ME  PLSQ  QUAOR AT  1 C  _ _ _ 

C  =  COEFFICIENTS  OF  5TH  ORDER  PLSQ  USED  TO  FIT  TI ME/TEMP .D ATA 
_  . RTENP  RECIPROCAL  .  ABSOLUTE  .TEMPERATURE _ ... 


_ EQ  U  I V  AjL.E  N.C.E.  (TLU  LjJTLOIT  1 )  l,i  ID , JZ I  1 )) , ( OATE1 , JZ  l  2  Lit  1 D ATE 2 , J Z (3 ) ) _ _ 

1,  (CCMWJZl  A  )  )  ,  (COM2, JZ(5)  )  ,  (COM3, JZI6)  )  ,  (COM4, JZ<7)  ) , < T INT , JZ < 8  )  > 

- DIJCEhSIOM _ ljjLSSHGI_LWIilJLQJ_LxTTIW  I IC  1 1  ,  DWDT  1 1 0 1 ) ,  Y  ( 1 20  )  ,C  I  6  )  .  A  ( I  ) , _ _ 

1W(55CC).,T(5  500),Z(42),X<36)  ,  PLOT  (  50 , 1 1  0  )  ,  TDER  (  1 01  )  ,  TPOLY  (  10  1  )  , 

_ LRT  E  M  P Ll.QJ.  J  J.E1  JJl  J  j  D_L6  I  j_J  Z  (.8  J _  _ _ _  _ 

INTEGER  DUMMY 

_ _ CATA  DUMMYZ.4JHZE.ro/.  _ _ __ _ _  _  _  _ _ 

98  CALL  READ(3,JZ,8,J)  "  '  2 

- IF  (  IC  ^EQ.L UMJLYJ.GJG  TJL  JJ.it 

CALL  WRITE(2»JZ»8) 


...9J  .CALL  _R..E  A  D  (  2  j  X  ,  18,J  )._ 
IF  (J-l)  111,112,113 


111  CALL  WRITE  (2»X»L8)  .  .  ._ 

13 

GO  TC  97 

n  ?  cam  r.i  csf  i  ?, ?i 

16 

GO  TC  98 

113  WRITE(6,11S  1JZI1) . . 

18 

115  FORMAT ( 2X, A5, 39HREDUNDANCY  ENCOUNTERED  -  RUN  TERMINATED) 

_ STEP  .  .  _  .  . .  ... 

11  A  READ! 5, ICCC ) ID*D A TE 1 ,DATE2 ,CGM1 , COM2 »CCM3, COM4 , TINT 
— I.CDC _ EJDiLMA  LIB  X.A5.1X«A6,A2.2X.3Aft.A3.GX,F6.4) 

19 

CALL  WRITE(2,JZ,8) 

21 

21 
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TGTAPE 

TC-TAPE 


09/17/7C 


-  EFN  SOURCE  STATEMENT  -  IFN(S)  - 


L  =  C _ _ _ 

ASSIGN  2C  TC  IEOF 

CALL  EOF  (  I  EC  F  )___ _ _ _ 

L=L  +  1 

N=18*L _ _ _ _ _ 

IF(N.GT.55CC)GC  TO  400 

M=N-1 7  _ ___ 

READ!  UT0C2  )  (  W(  I  )  ,T!I)  ,I=;M,N) 

FORMAT  (3fc(F6.2»lX)  I _ _ 

I  CALl  V»R  ITE  {  2  »W(  M  J  ,  18  ) 

C ALL  W R1TE(2,T(M),  18) _ _ 

'  GO  TC  10 

L=L— 1 _ 

N=18*L 

CALL  CLOSE  <_2  »_2J _ _ _ _ _ 

CALL  WRITE ( 2  » DUMMY  ,8  ) 

JLK  =  L-5_ _ _  _ _ _ 

CO  38 C  K= 1 ,  LK  ,  5 

WRITE  (  6.4)  IC.K _ ___________ _ 

"FORMAT (lHL, 10X  ,A5 , 5X »6HREC  NO , 1 3 2/38X, 3 1 1HW , 

_ KT =  K  +  4 _ _  _  _  _ 

CO  370  J^K »  K  T 

_KB^_18*J-17 _  _ _ 

KE  =  18*J— 12 
DO  365  I  =  KB » KE 


i  WRITE(6,5)W( I) ,T(I),W( 1  +  6  ),T(I+6),W( I  +  12),T< 

F0RMAT(35X,61F7.2,4X1) _ _ 

WR  ITE (6,381  ) 

FORMAT!  1H0 ) 


t  CONTINUE 

:  CONTINUE _ _ _ 

LR  =  L-U/5)*5 

IF ! LR . EC . 0 ) LR-5 _ _  _ 

M=L-L  R  + 1 

WR I TE I  6,4 )  I C  ,M _ _ _ _ _ _ 

CO  385  J  =  M,L 

KB=1P*J— 17 _ __ 

KE=16*J-12 

„£1J  _3S5_J3K_E_,  KE _ _ _ _ 

WR  I TE ( 6  *  5 )  W (I),T( I) ,W( 1+6) ,T( I+fi),W( 1  +  12)  ,T 

WRITE  I  6,381  ) _ _ 

CONTINUE 

. vl  J = ..« C  L  C  f.E  UL6X1  Ml _ 

LL  =  MAXO ( J  J ,  1 C  ) 


11X.1HT,  11X) ) 


(  1  +  12  ) 


JJ  =  1  PERCENT  OF  NO.  OF  DATA  SETS  READ  IN 
LL  =  NO.  OF  CURVE  F_IT_  PO  INTS _  J  L  4JJR__=_NN  ) 


WRITE  (6.3CC0) _ _ 

SCOC  FORMAT  ( 1F1 ) 

WRITE  (6,3  C  5  0 )  I_D ,  CATE  1,  D  ATE2  ,C0M1  .COM2, COM  3,  C0M4_ 
)050  FORMAT  (4X,  A5,8X,A~6,A2,  lCX,A6t.A6,A6,  A3//I 

„ _ WRITE!  fe.306C)TINT _ _ _ _ _ 

i C60  FORMAT ( 10X.14HTIME  INTERVAL^ F6.4 ) 


3 17C  FORMAT  (10X.25HNC  CF  PTS  IN  CURVE  FIT  =  ,12) 
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TCTAPE  09/17/7C 

TGTAPE  -  EFN  SOURCE  STATEMENT  -  IFN(S)  -  .  .. 

WRITE  (A, 30)0)  N 

106 

3CIC 

FORMAT  (10*21  ETOTAL  NO  CF  POINTS  =  ,14) 

c 

Sc 

K  =  POLYNOMIAL  ORDER .NEEDED  FOR  PLSQ  SUBROUTINE.  LIST  =  0  FOR  NO  ERROR 

c 

f. 

ANALYSIS  CF  PLSQ 

D  =  T0TA1  WEIGHT  LOSS 

c 

LINDA  =  1  .  _  _ _ _  _ .  ....  ...... 

D  =  Ml)  -  WIN) 

B  (  1  )=-6.8S5309E-6 

B( 2  )  =  3.5219C5E—  4 

B( 3 1 -  —  7 . 78  2  8  05F—  3 

B  (  4  )  =  9.75327E—  2 

81 5)=-7.656367E-l  _  _ _  _ _  _ _  . 

B(6)=3.943215EC 

B(7)=-1.367422E1 

B  ( E)=3.2741fi6El 

R ( 9  )  =—  5.749016F1  _ 

B (  10)=1.819171E2 

fill 1)=3.RI 2777F-?  _ 

DO  55  1=1, N 

urn  =  ioo.-iico.*(wm-MN))/D)  ...  . . . . 

TI(I)=(2.*FL0AT(I)-l.)*TINT/60. 

POL  Y=  Bill  _  _  .....  _ 

3  or. 

CO  3CC  J-2,11 

POL Y=POLY*T ( I )/10.+B( J) 

55 

T ( I  )  =  PCLY 

CONTINUE 

C 

C 

CURVE  FIT  OF  TIMF  AND  TEMPERATURE  DATA 

c 

K  =  5 

LIST  =  0 

CALL  PLSQ(TI.T.N.K.D.LIST,EMAX,ERMS,EMEG> 

129 

5100 

WRITE  16*5100)  EMAX 

FORMAT  (10X.17HMAX  TEMP  ERROR  =  .F10.6) 

130 

_ 5200 

WRITE  (6.52C0)  ERMS 

FORMAT  ( 10X  *  30HTEMP  ROOT  MEAN  SQUARE  ERROR  =  .F10.6) 

131 

5300 

WRITE  (6,5300) 

FORMAT  (10X.15HTEMP  POLY  COEFF) 

132 

5400 

WRITE  16,5400)  (D(I),I=1,6) 

FORMAT ( 13X.F1 2.6) 

133 

c 

c 

START  MAJOR  LOOP 

c 

CO  ICO  NW  =  1.99 

58 

II  =  LINDA-1 

WW(NW)  =  FI OAT(NW) 

C 

c 

SCAN  WF  TGHT  DATA  FOR  ONE  CLOSE  TO  BUT  JUST  GREATER  THAN  ONE  PERCENT  WEIGHT 

c 

s 

LOS  St.  II  =  INDEX  OF  THAT  POINT 

CO  6C  I  =  L  INEA.N 

— _  II  =  LL*1 _ 

IF  (W(  l  1.GT.WW1NW)  )  GO  FC  70 
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TGTAPE _  - _ EFN  SOURCE  STATEMENT  -  IFN(S)  - 


6C 

CONTINUE 

c 

7  C 

LINDA  -  II-ILL/2) 

c 

c 

LINDA  =  INDEX  CF  FIRST  DATA  TO  BE  USED  BY  PLSQ 

DO  8C  J=1,LL 

JI  =  L  INDA+J-1 

T I ( J)=! 2. AFLOAT ( J I )-2. ) ♦TINT/ 60. 

Y ( J )  =  W(JI) 

C 

8C 

CONTINUE 

C 

C 

CURVE  FIT  CF  TIME  AND  WEIGHT  DATA 

K  =  2 

LIST  =  C 

CALL  PLSQU I, Y,NN,K,C,LIST,EMAX,ERMS,EMEG) 

KK  =  1 

168 

C 

C 

START  LCOP  TO  CHECK  FOR  BAD  INPUT  DATA 

C 

DO  81  J=  1  *  L  L 

C 

JI  =  L  INDA+J-1 

c 

c 

WE  =  WEIGHT  CALCULATED  FROM  POLYNOMIAL 

c 

WE  =  Cf  1)*TII  J  )**2+X(2)*TI  (  JI+CI3) 

c 

c 

COMPARE  CALCULATED  AND  ORIGINAL  DATA 

IF  (ARSIWE-WIJl ] ).GT .1.)  GO  TO  82 

GO  TC  81  _  .  _  ..  _  ..  .  . 

82 

4000 

WRITE  1 6 » 4  C  CO )  JI,W(JI),WE 

FORMAT  ( 1  OX » 9HA  T  PT  NO  ,I4,10H  WEIGHT  =  ,F5.1,13H  REPLACED  BY  , 

182 

c 

•F5.1) 

c 

X 

REPLACE  BAD  DATA  BY  CALCULATED  VALUES 

W(JI)  =5  WE 

KK  =  2 

81 

CONTINUE 

GO  TC  (83.S81.KK 

c 

JL. 

CHECK  FOR  IMAGINARY  ROOTS  IN  SOLUTION  OF  QUADRATIC 

c 

83 

SCREW  =?  CI2)*CI2)-4.0*Cm*ICC3)-WW(NW)} 

c 

IF  (SCREW. LT. 0. 0 >  GO  TO  90 

c 

x_ 

USE  REAL  RCCT  TO  DETERMINE  TIME  CORRESPONDING  TO  EACH  PERCENT  WEIGHT  LOSS 

c 

TNW(NW  )  =  ( SQRTICI 2 J*C(2)-4.0*C( 1 ) * (C ( 3 > -WW l NW ) ) ) -C ( 2 ) ) / ( 2 ,0*C ( 1) ) 

c 

0 

CWCT  =  RATE  OF  WEIGHT  LOSS 

194 

DWCT(NW)  =  2.0*C  (I)*T.NW(NW)  M  C(2) 

GO  TC  100 

90 

TN W(N W  1  =  12. ♦FLOAT)  II )-4 . ) *T INT/6C. 
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_ _ TOT  AP  E 

nwrT(MW)  =  c.o 


09/17/7C 


_ EFN _ SOURCE  STATEMENT  -  IFN(S) . 


C 

c 

c 


HR  1 16  pLT  I  CENT  IFICAT1GN  ANO  IQCATI  ON_OF._  BAJ).  DA  TA _ _ _ 

_WR  I T E  <  3 U  C  )  NW 1 1 1 ,  TNW  ( N W  )  ,.W(  ID _ _ _ _ _ _ _ _  _ _ _ 

3160  FORMA  K2X,  17HSCREW  LESS  THAN  0 , 1CX, 3HNW=, I  3 , 10X ,3HI I  =  , 1 4 » l OX, 


100  CONTINUE 

_ WRITE  16.3110) _ _ _ _ _ _ _ _ _ . - 

3110  FORMAT (//3X» 11HWEIGHT  L0SS.6X  , 8HDWDT (NW > , 14X, 4HTEKP ,6X , 

_ *it_H  T £1  E.R  »JLLX j  SHKlf.M P.,.16X  , 4KTIME). _ _ 

CO  120  NW= 1,99 

_ CI=.IN.WXtiWJ - _ - - - __ - — - - - - — 


202 


.207 


TSTCR  -  D(  1 ) 

C  LOOP  TO  EVALUATE  TEMPERATURE  POLYNOMIAL  FOR  EACH  VALUE  OF  CT 

DO  2C0  D2 , 6 

_ 2.0  C  TSTCR.  5  TS1CR*CHID(I ).  - - - - - 

TPCLYINW)  =  TSTOR 

_ . TSTCR  ,=.  5«*D  1 1 ) . . . . . . . . . . . . - . . - 

DO  250  1=2,5 

250  TSTOR  =;  TSTCR*CT+TL0’aTIJ  j*D(  I) 

C  TDER  =  TEMPERATURE  DERIVATIVE 

_C _ RTEMP  =  RECIPROCAL  OF  ABSOLUTE  TEMPERATURE _ _ _ 

C 

. IDER.LNW)  .=  TSTOR . . . . . .  . . . . . . .  . . 

RTEMP  I  NW  )  =  l.C/(TPOLY(NW)+273.16) 

_ WRITE  (6,3120)  NW ,  DWDT  ( N  W  )  ,  TPOL  Y  (NH  )  ,  TDER  ( NW^RlEMPlNMlxO _ 

3120  FORMAT  t 6X , I  3 , 10X , El 2 . 5 , 7 X , F9 . 3 , 2E 1 5.5 , 5 X, F7. 2 ) 

.120  CCMINUE  _ _ _ _  _ _ _ _ _ _ _ - . 

STCER  =;  O.C 

C  CALCULATE  AVERAGE  TEMPERATURE  DERIVATIVE  (AVE) 

X _ 


2.26. 


CO  125  1=1,99 

_ _S.TJC.ER _=,  STCER  4  TDER  ( 1  ) _ _ _ 

125  CONTINUE 

AVE  =  STDEP  i 99,0. _ _ _ _ _ 

WRITE  16,3125)  AVE 


241 


,  E 1 5, 5J_ 


C 

X 

C 


SET  LP  DUMMY  POINTS  FOR  GRAPH  PLOTTING  SUBROUTINE  (GP) _ 


WW(  ICC)  =  X_. Q _ 

CWCT(ICO)  =  C.C 
_INW  ( ICC  ).^_INW.l  99.1. 


TDER (ICO)  =  O.C 
W«.LIC11_=.  1C  0  .C. 


CWCT(lCl)  =  C.C 
.  _T  Mte.LLQ  L) ...  =_.IN  HL9.91.4_ 
TDERI101)  =  TDERI99) 
WRITE.  IX, 3.C.C0J _ 


ZhZ 


WRITE  (6,3130 


ID 


243 
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TGTAPE  -  EFN  SOURCE  STATEMENT  -  IFN(S) _ - 


3  13C 

FORMAT  (  1CX  ,  19FDWDT  VS  WEIGHT  L0SS.20X.A5) 

L  =  3 

LS  =  5 

LW  =  1C1 

LN  =  50 

y  =  ICE 

DATA  A/EH./ 

C 

JN  =  1 

C 

C 

PLOT  GRAPH  GF  RATE  OF  WEIGHT  LOSS  AGAINST  PERCENT  WEIGHT  LOSS 

CALL  GP  (WV.DWCT.L.LS.M, JN , LW , LN , A,  PLOT ) 

WRITE  i  6. 3CC0 ) 

250 
25  1 

3  1 A  0 

WRITE  (6,3140)  ID 

FORMAT  ( 10X  ,  19FWEI GHT  LOSS  VS  TIME.20X.A5) 

252 

C 

c 

c ' 

_PLCJ_ GRAPH  CF  PERCENT _ WEIGHT  LOSS  AG A  INS  T  TIME 

CALL  GP  <TNW,WW,L,LS»M,JN,LW.LN,A,PLOT) 

253 

WRITE  I6.3CCC) 

WRITE  16,3150)  ID 

254 

255 

3150 

C 

FORMAT  (  1CX ,  l  2HTDER  VS  TIME.20X.A5) 

C 

C 

PLCT  GRAPH  CF  TEMPERATURE  DERIVATIVE  AGAINST  TIME 

CALL  GP  ITNW.TDER.L.LS.M.JN.LW.LN.A.PLCT) 

GO  TC  500 

256 

400 

6000 

WRI TE(6,6CCC ) 

FORMAT! 10X.48HNUMBER  OF  DATA  POINTS  EXCEEDS  NUMBER  DIMENSIONED) 

258 

50C 

STEP 

PROGRAM  2 


This  program  provides  print  out  of  rates  also  rate  and 
temperature  on  punched  cards  for  use  in  Program  3. 
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TGTAPE _ 2 _ EJFN _ SOURCE  STATEMENT  -  IFN(S)  -  _ 


c 

c 

PROGRAM  2  CUT  OFF  AND  RATE  DATA  FOR  A  SERIES  OF  RUNS 

PROGRAM  TO  RE AC  A  SERIES  OF  RUNS  FROM  THE  MASTER  FILE  ,  APPLY  THE 

c 

c 

APPROPRIATE  CUT  OFF  VALUE,  AND  OUTPUT  THE  RATE  OF  WEIGHT  LCSS  AT 

1  PER  CENT  INTERVALS  ON  CARDS  FOR  USE  IN  THE  ARRHENIUS  PROGRAM 

c 

c 

INPUT  TEMPERATURES  ARE  FITTED  TO  A  FIFTH  DEGREE  POLYNOMIAL  USING  A 

LEAST  SQUARES  SUBROUTINE  (PLSQ). 

c 

c 

WEIGHTS  CORRESPONDING  TO  SHORT  TEMPERATURE  RANGES  ARE  FITTED  TC  A 

CUACPATIC  BY  PLSC. 

c 

c 

INPUT  WEIGHTS  DIFFERING  FROM  FITTED  LINE  BY  MORE  THAN  ONE  PERCENT  CF  THE 

TCTAL  WEIGHT  LOSS  ARE  REPLACED  BY  THE  CURVE  FIT  VALUE, 

c 

c 

CUT  PUT  DATA  IS  PUNCHED  ON  TO  CARDS  FOR  FURTHER  PROCESSING  (TO  CALCULATE 
ACTIVATION  ENERGY  ETC). 

c 

c 

W  =  WEIGHT  DATA  POINT  READ  OFF  TAPE 

T  =  TEMPERATURE  DATA  POINT  READ  OF  TAPE ( IN  MV.) 

c 

c 

c 

c 

— 

T I  =T IME  DATA  POINT  CALCULATED  FROM  TIME  INTERVAL  AND  NO  OF  DATA  POINTS 
TINT=TI ME  INTERVAL 

WW  -WEIGHT  LOSS  AT  1  PER  CENT  INTERVALS 

CWCT  =  RATE  OF  WEIGHT  LOSS  AT  1  PER  CENT  INTERVALS 

— 

c 

c 

TGER  =  HEATING  RATE  AT  1  PER  CE'NT  INTERVALS 

TPDLY=TEMPERATURE  CORRESPONDING  TO  EACH  PER  CENT  WEIGHT  LOSS , CALCULATED  BY 

c 

c 

PLSC 

TNW=  TIME  CORRESPONDING  TO  EACH  PER  CENT  WEIGHT  LOSS 

c 

c 

PLOT  =D  IMENS IONS  FOR  GRAPH  PLOT  SUBROUTINE 

6  =  COEFFICIENTS  OF  10TH  ORDER  POLYNOMIAL  FITTING  TEMP/EMF  DATA  FOR 

c 

c 

PLATINUM*PLATINUM  10PER  CENT  RHODIUM 

C  =  COEFFICIENTS  OF  WEIGHT/TIME  PLSQ  QUADRATIC 

c 

c 

C  =  CCE FFIC I ENTS  OF  5TH  ORDER  PLSQ  USED  TO  FIT  T I ME/TEMP. D ATA 

RTEMP  =  REC  I PROC AL  ABSOLUTE  TEMPERATURE 

c 

ECU  I  VALENCE  IT  I  (U,  PLOT  (1)  ) ,  (  I D,  JZ  ( 1 )  ) ,  ( DATE  1 ,  J  Z  ( 2  ) )  ,  (DATE2  ,  JZ  (3  ) ) 

1,  (CCM1, JZ(4) ), (COM2, JZ(5)), (COM3, JZ<6)), (COM4, JZ (7)> , (TINT,JZ(8) ) 

DIMENSION  TI(5500),WW< 101),TNW( 10 1 ) , DWDT ( 101 ) , Y ( 12C ) ,C(6) ,A(1) , 

1WI55C0) »T(5500),Z(42)»X(36)»  PLOT ( 50, 1 1G I ,TDER ( 101 ) ,TP0LY<101) , 
1RTEMF(101),B(11),D(6),JZ(8) 

DIMENSION  I C  A ( 12 ) »  OCA ( 12  ) 

INTEGER  DUMMY 

CAT  A  CUMMY/4FZER0/ 

N  F=  1 

15 

1010 

REAC(5,  1010) IDA (NF),CCA(NF) 

FCRMAT(8X,A5,1X,F5.3) 

2 

IF(ICA(NF).EC.CUMMY)GO  TO  20 

NF=N  F  +  l 

20 

GC  TC  15 

NF=NF-1 

99 

CO  5  IN F= 1 ,NF 

CALL  READ(3,JZ»8,J) 

16 

IF( IC.EC.CUMMYJGO  TO  101 

IF( IC.EQ.ICAI  INFJJGO  TO  97 

98 

C  ALL  RE AD(3»W» 18 ,  J  ) 

I F ( J-l ) 98 ,99 » 98 

26 

101 

4500 

WR IT  E ( 6 ,4500 ) 

FORMAT) 10X.25HSEARCH  EXCEEDS  VALID  FILE) 

29 

97 

STOP 

CC=CCA ( INF) 

L=Q 
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25 

M= 18  *L+  1 

CALL  READ(3,W(M) ,18; J) 

IF { J- 1 » 10,51,30 

34 

10 

CALL  READ(3,T(N),18, J) 

1 F  (  J- 1 )  40, 50, 30 

39 

“  50 
1000 

WRITET6,1DW)XC,L 

FCRMAT(2X,21H0UT  OF  PHASE  DATA  IN  ,A6,18H,  AFTER  RECORD 

NO  ,13) 

43 

40 

TTOF 

L  =  L+  1 

N=18*L 

GC  TC  25 

30 

6 

WRITE  16 , 6)TZ( I ) ,1=1 , 42 ) 

FCRMATl  1H1, 10X.39HTAPE  READ  ERROR  IN  THE  FOLLOWING  RECORD/ {8X , 14 ( A 

47 

16 , 2  X  )  > ) 

STCP 

^1 

J J= .0 iO ‘FLOAT  ( N ) 

LL  =  MAXO ( J J , 10  ) 

C 

c 

JJ  =  1  PERCENT  OF  NO.  OF  DATA  SETS  READ  IN 

t 

c 

LL  =  NO.  OF  rORVFTTT POINTS  (LATER  =  NNT~ 

3000 

WRITE  (6,3000) 

FORMAT  (1H1) 

54 

3050 

WRITE  (6,3050)  1 0 ,HffTFT,0£TE2, COM  1, COM2, COM3, COM4 
FCRMAT(4X, A5,8X, A6,A2, 10X, A6,A6,A6,A3//> 

55 

3060 

WRITE(6,30  6CT)  TINT 

FGRMATt 10X, 14HTIME  INTERVAL3 , F6 . 4 ) 

56 

3170 

WRITE  (6,3170)  LL 

FORMAT  ( 10X , 25HN0  OF  PTS  IN  CURVE  FIT  =  ,12) 

57 

3010 

WRITE  (6,3010)  N 

FORMAT  ( 10X , 21HT0TAL  NO  OF  POINTS  =  ,14) 

58 

NN  =  LL 

WR IT  E ( 6 , 3020  )  DC 

59 

3020 

C 

FCRMATl 10X.9HCUT  OFF  =,F5.3) 

C 

c 

K  =  POLYNOMIAL  ORDER, NEEDED  FOR  PLSQ  SUBROUTINE.  LIST  = 
ANALYSIS  OF  PLSQ 

0  FOR  NO 

ERROR 

c 

c 

C  =  TOTAL  WEIGHT  LOSS 

UNCO  =  1 

C  =  Ml)  -  WIN) 

Etl)=-6. 885309 E- 6 

E(2)=3.52l905E-4 

e<3)=-7.783805E-3 

E(4)=9.75327E-2 

E(5)=-7.656367E-1 

E ( 6  )  =  3. 9432 15  EO 

E(7)=-l  .367422E1 

E ( 8 ) =3. 274 186E 1 

Et9)=-5.749016E1 

B ( 10 ) *  1 • 8 19 17 1E2 

etll)=3.812777E-2 

CC  55  1  =  1, N 

MI)  =  100.-(100.*(W(  I)-W(N)  )/0) 

W ( I ) =W (  I)/OC 

Tit  I  )  =  ( 2. * FLOAT ( I )  —  1 . ) *T INT /60. 
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FCLY  =  e(_l)_ 
cc  SCO  J=2, ii 

300  PCLY=POLY*T( I)/10.+B(J) 

T (  I  )  =  POLY 
55  CONTINUE 
C 

C  CURVE  FIT  OF  TIME  AND  TEMPERATURE  DATA 

c . .  '  . . 

K  =  5  _ 

LIST  =  0 

CALL  PLSCITI ,T,N,K, O.LIST, EM AX , SRMS , EMEC ) 

WRITE  (6,5100)  EMAX 

5100  FORMAT  (10X.17HMAX  TEMP  ERROR  =  ,F10.6) 

WRITE  (6,5200)  ERMS  " 

5200  FORMAT  ( 10X,  30FTEMP  ROOT  MEAN  SQUARE  ERROR  =  , F 10 . 6 ) 

WRITE  (6,5200) 

5300  FORMAT  (10X.15HT5MP  POLY  COEFF) 

WRITE  (6,5400)  (0(1),  1  =  1, 6) 

5400  FCRMATl  13X, E12.6  )  _ 

C 

C  START  major  LOOP 
c 

00  100  NW  =  1,99 
58  II  =  LINCA-1 

WW(NW)  =  FLOAT ( NW ) 

C 

C  SCAN  WEIGFT  CATA  FOR  ONE  CLOSE  TO  BUT  JUST  GREATER  THAN  ONE  PERCENT  WEIGHT 

C  LOSS.  II  =  I N 0 E X  OF  THAT  POINT 
C 

CC  60  I  =  L  I N C  A  ,  N 
II  =  II  +  l 

IF  (W(I  )  .  C-T  .  WW  ( NW  )  )  00  TO  70 
60  CONTINUE 
70  LINOA  =  I  I-( LL/2 ) 

C 

C  LINOA  =  INCEX  OF  FIRST  OATA  TO  BE  USED  BY  PLSQ 

C  . . 

CC  80  J  =  1 ,  LL 
Jt  =  LINCA+J-1 

T  I  ( J)  =  ( 2.  *  FLO AT ( JI )-2. )*TINT/6C. 

Y  (  J  )  =  W  (  J  I) 

80  CONTINUE 
C 

C  CURVE  FIT  OF  TIMS  AND  WEIGHT  DATA 

C 

K  =  2 
LIST  =  0 

CALL  PLSQITI ,Y,NN,K,C,LIST,2MAX,ERMS,EMEQ) 

KK  =  1 
C 

C  START  LOOP  TO  CHECK  FOR  BAD  INPUT  DATA 

C 

CC  81  J  =  1 , LL 
JI  =  LINCA+J-1 

C  _  _ _  _ 

0  WE  =  WEIGHT  CALCULATED  FROM  POLYNOMIAL 
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CTTT*t  I  <  J  ) **2-»C 1 2 ) »TI  ( Jl'-fCTJT  . ~~  . . 

CCM FA  RE  CALCULATED  AND  ORIGINAL  DATA . . 

IF  (ABS(WE-W(Jl)).GT.l.)  GO  TO  82 
CC  TC  81 

82  WRITE  ( 6 * AOCO~ »  jl,W(JI),WE"  ''  “  .  .  .  . 132 

00  FORMAT  (10X.9HAT  PT  NO  , 14,  1CH  WEIGHT  =  ,F5.1,13H  REPLACED  BY  , 

•  F  5  .  1  ) 

REPLACE  BAD  DATA  BY  CALCULATED  VALUES 

WUTf  =  ~WE~ .  . . ""  . . .  ’  . . 

KK  =  2 

81  CONTINUE  . 

CC  TC  (83,58), KK 

CHECK  FOR  IMAGINARY  ROOTS  IN  SOLUTION  Of  QUADRATIC 

83  SCREW  =  C(2)*C(2)-4.0*C(1)*(C(3)-WW(NW)) 

'  IF  (SCREW.LT. 0.0)  GO  TO  90 


US:  REAL  ROOT  TO  DETERMINE  TIME  CORRESPONDING  TO  EACH  PERCENT  WEIGHT  LOSS 
"  TKwTNWf  =  LS  grTTC  ( 2 ) *C ( 2 T-47c*CTiT* ( C(3  ) -WW ( NVOTI-C  ( 2 )T /T2 .0*C ( IT )_ 


CWCT  =  R  AT  EOF  WEIGHT  LOSS 

144 

CWCT  (NW  )  =  2.0«C(  i)Vf'NW'(NW)  +  C(  2) 

GO  TC  100 

90  TNW  ( NW )  =  (  2.*  Fid  AT  (  f  I  )-4  ,)  #tTKt /60. .  .  . . 

CWCT(NW)  =  0.0 


WRITE  OUT  IDENTIFICATION  AND  LOCATION  OF  BAD  DATA 


__ WRITE <6 , 3160)NW, II,TNW(NW)  ,W(  II) 

160  FORMAT  ( 2X , 1 7HSCREW'  Tf SS  THAN  C, 10X , 3HN W= , I j , 10X , 3H 1 1  =  , 1 4 , 1 CX , 
. 2HT=,F6.2,10X,2HW=,F5.1) 

100  CONTINUE 

WRITE  (5,3110) 

110  FORMAT ( //3X , 1  lFWE  IGHT  LOSS , 6 X , 8HDW0T ( NW ) , 14X ,4HTEM P,  6 X , 

•  4HTCER, 11X,EHRTEMP,  16X.4HTIME  ) 

CC  120  NW  =  1 , 99  . . . . .  '  . . "  “  . . 

CT  =  T  NW ( NW  ) 

t  stop  =  cm 


152 


156 


LOOP  TO  EVALUATE  TEMPERATURE  POLYNOMIAL  FOR  EACH  VALUE  OF  CT 

CC  2C0  I  =  2,6  . .  . “  '  . . . . 

200  TSTCP  =  TST0P*CT+0(  I  ) 

TPCLY(NW)  =  TSTOR 
TSTCR  =  5.*D( 1  ) 

CC  250  1=2,5 
J  =  6-1 

250  TSTCR  =  TST0R*CT  +  FL0 AT ( j  )*0 (  I  ) 
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.  f C £ R  =  TEMPERATURE  DERIVATIVE  -  -  -  .  -  -  -  . 

RTEMF  =  RECIPROCAL  OF  ABSOLUTE  TEMPERATURE 

TQER(KW)  =  TSTQR 

RTSMP(NW)  =  1  .0/ ( TPOL Y ( NW  )  +  27  3 . 1 6  ) 

W R I TF  (6,3120)  NW .DWCT (NW), TPOLY(NW) , TOER (  N  W )  ,  RJ  E MP(NW) ,CT _  _  _  _  174 

120  FCRMAT  ( 6X , I 3 , 10X, E l 2, 5 , 7X , F9 .3, 2E 15 . 5 , 5X , F7.2 ) 

120  CCMINUE  _ 

STEER  =  O.n 

CALCULATE  AVERAGE  TEMPERATURE  DERIVATIVE  (AVE)  .  . 

*1  =  j.7gq'  ’  ~  _ 

_STCER  =  STDER  4  TDER(I) 

125  CCNTINLIE  " 

AVE  =  STDER/99.0  _ _ _ _ 

WRITS  (  6,3125  )  AVE  .  ~  ~  '  186 

125  .  FCRMAT  t//10X,27H  AVFRAGE  TEMP  DERIVATIVE  =  .E15.5)  _ _ _ -  . . 

SET  UP  DUMMY  POINTS  FOR  GRAPE  PLOTTING  SUBROUTINE  (GP) . .  .. 

WWt  IOC)  =  0.0 
CWCTUOO)  =  b.c 

JNWllCOJ  =  TNW(_99)  __  _  _  _  _  _  _  _  _  _  _ 

TCEPt  10  0 )  =0.6 
V.  W  C  101)  =  100. U 
CWCT(lOl)  =  0.6 
TNW(lOl)  =  TNW ( 99 ) 

TCEP(lOl)  =  TCERI99) 

WRJTE  (6.3000) _  .  _  _  _  _  187 

WRITS  (6,3130)  ID  188 

13C  FCRMAT  ( 10X . 19HDWDT  VS  WEIGHT  LOSS.20X.A5) 

L  =  3 

LS  =  5  ___  _  _ 

LW  =  10  1 

LN  =__50 _  __ _ _  __  _  _  _ 

M  =  101 

CATA  A/ 1 H . /  _  _  _ 

JN  =  1 

FLCT  GRAPH  OF  RA T  E  0  F ' W E I G H T  LOSS  AGAINST  PERCENT  WEIGHT  LCSS 

WU'^TWW,  CWCT.L.LS.M,  jT.ThtLN'iA/Plofi  . . .  .  195 

WRITE  ( 6 » 3000 )  _  _ _  __  _  _  196 

CWCT(IOO)  =  C.C 

_ T PC L Y  1100)  =.  0.0 . . . . . .  . . . . . . . . . . .  . 

_ PUNCH  OUTPUT  CARDS  CONTAINING  PERCENT  WT.  LOSS(NW)  THEN  THREE  PAIRS  CF _ _ 

TEMPERATURE  ANC  RATE  0 F  W E I G F T  LOSS  DATA 

CC  15  ON  W  =  1 »  Too » 3  ”  '  """  . 

PUNCH  5000. IC.NW. DWCT(NW ) , TPOLY ( NWJ ,DWDT (NW4l ) , TPOLY (NW+1J  »  _ 

*  CWCT  (NW4  2  ) ,  T  POLY  ( NW  +  2) . . . .  ”  '  198 

5COO  FCRMAT  1 1 X ,  A5 , 1 A .  E 13 .5 ,  F6. 1 ,  E  13. 5.  F6 . 1 ,  E 13. 5  ,F6 . 1 ) _ _ _ _ _ _ _ 

ISO  CONTINUE 

TGTAPc .  .  .  . 09/17/70 

TGTAPE  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 


5  CCNT1NUE 
'  STCF 
ENC 
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PROGRAMME  TC  DETERMINE  ./TGA  PARAMETERS  BY  FRIEDMANS  METHOD _ _____ _ _ _ _ 

PROGRAMME  ACCEPTS  DATA  CARDS  HAVING  THREE  SETS  C3F  DATA  PER  CARD. 

LAST  CARD  Of.  EACH  DECK.  MUST  .HAVE  A  ONE  IN  COLUMN! .  LAST  CA_RD  GF  _  ___  _ _ 

LAST  DECK  FOR  CNE  POLYMER  SYSTEM  MIST  HAVE  A  TWC  IN  COLUMN  1  INSTEAD 

TO  RUN  A  SECCNC  SET  CF  DECKS.  PUNCH  A  CARD_WlTH  A.  THREE  INJSCLUMN  1  _ 

AND  PLACE  BETWEEN  SETS 

_ AT  THE  END  OF  ALL  CECKS  PLACE  A  BLANK  CARD  THEN  AN _ 4E OF  ______ _ _ _ 

SYMBOLS  __DWBT  =  RATE.  C_F  WEIGHT  LOSS*  RTEMP  »  RECIPROCAL  0_F_AB  SOLUTE 
TEMPERATURE,  PATE  =  LCG  RATE  OF  WEIGHT  LOSS*  SLCFE  =  SLOPE  OF  ARRHENIUS 
PLOT,  P  R£  X  tP^RE  - EXPONENTIAL  FACTOR,  PLOT  =  DIMENSION  OF  GP  SUBROUTINE  _ 

ACTE  =  ACTIVATION  ENERGY,  X  AND  Y  REPRESENT  DATA  TREATED  BY  GP 

JPQLY  =  _  I NPU I  _T  E  M  Pi  R  A I U  P  E  S ,  ID  »  I C  ENT  IFICATIDN  ,  A  =  NOi.  OF  SYMBOLS  IN  GP _ 

AA  =  PERCENT  WEIGHT  LOSS,  AFW  =  FUNCTION  FROM  FRIEDMANS  EQUATION 

FW  =  AVERAGE  AFW,  88  *  LOG  (PERCENT  RESIDUE),  WF  «  AVERAGE  AFW _ _ 

DI MENS  I CN  DWDT ( 100,10), RTEMP ( 100, 1C  ), RATE ( 100, 1 C  ), SLOPE (100 ) , 

.PREX(IOO) ,PLCT ( 50, 100 ) , ACT E ( 1 00 ) , X ( 10),Y( 10), 

_  ..TPOL.Y  < JL  C_C  ,10)  ,  LDUO  ) ,  A  ( 1  )j  A  A  <  101)*A  FWf  10 ) .  EW(  10  C  ) ,  B8  (  95  ) ,  WE  (95  )  , _ _ _ 

•  SPS(IOO) , SOS (100) , SO  I  (  ICO) ,B(8) 

1  READ  (5,1000)  I  GfCOMl ,  CCM2 ,  COM3 » COM 4, COM 5 ,  C0M6  ,_CCM7.f  COMB  .  . . _ . 1. 

WRITE  (6,3000)  3 

W R ITE  ( 6 , 1100  )  I Gi.CQMl , COM2  ,C0M3, COM4,  COM  5,C0M_6  ,CQM7, CO M.8  _ _ _ 4 

2  J  =  0 

_  1Q.J . =_.JU _ 

_ START  LCCP  TC  READ  IN  DATA _ _ _ 

...  00  20  .Nh  *  _1_,  91,  3 _ _ _ 

LB-J_  =1.  IN.  .CCLUMN.1  CF  LAST  CARD  Of  A  CECK,  LAST  CARD  OF  LAST  DECK  FOR  _ „ 

ONE  POLYMER  SYSTEM  NEEDS  LBJ  =  2. 

READ  (5,1200)  L B J , I D ( J ) , I W , CW DT ( NW , J ) , TPCL Y (NW , J  ) , DWDT ( NW  + 1 , J ) , 

•  TP0LY.1NW+JL,  J_1  ,CWDT(NW»2,J),TP0LY(NWY2,J) _ _ _ _ 

_ _CH ECK  THAT..  JN EUt_.CARCS  _A_R.E__IN.XQN SECUTIVE  ORDER  .  . 

10 

_  IE  _(  I W-NW ! 3  *4  ,3 _ _ _ _ _ _ _ _ _ 

3  WRITE  (6,1900)  NW,IO(J),IW  21 

4  AA(NW)  =  FLCAT(NW) 

_ .AA  LN.W+ 1  )__=_EL  CAIjLMHIU _ _ _ _ _ . _ 

AAINW+2)  =  FLC  AT  (  NW-K2  ) 

_ IF .  JLBJ.EQ.l)  GO  JH  .10 _ _ _ _ _ _ _ _ _ 

20  IF  (LBJ.EQ.2)  GO  TO  25 

. .  WRITE  .  LIST  GF  .RUN  IDS _ _ _ _ _ _ 

.  _  WRITE  (6j1800MI_D(IJl,_1*1.J) _ _ _ _  _ _ _ _ _ _ _ _ _ _ 

.  CHECK  FGP  AT  _Lf  AST  _T_hE  EE  DATA  OECKS _ _ _ _ _ _ _ 

36 

- IF  (J-3).  .30,35,25 _ _ _ — - - 
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30  WRITE  (6,2000)  44 

_ .t  -G B TC  3 QC . __ __ . — _ _ _ _ _ _ _ _ _ _ _ _ 

35  WRITE  16,1500)  46 

. ...  _  -T5UM  .  =  .Q _ _ _ _ _ 

N  =  0 

SPREX  .a  . 0.0 _ _ _ _ _ _ _ _ _ 

_ START  -LHC-E  -HI ..CALCULATE -LEASt-^SQUARfS^LlNE-nF-LCCl&ATEJ...  _ _ _ 

_ DO  45  NW  =  4,58  _ _ _ _ _ _ _ _ _ . _ ___ _ — . . . . 

SUMXX  =  C 

. .  SUMY.Y  =  .a _ _ _ _ _ , _ - _ _ _ _ _ _ _ _ _ - 

SUMX  =  0 

__ S UMY._n._a _ _ _ _ _ _ _ _ _ ___ _ _ _ _ __ _ _ _ _ _ ___ 

SUMX Y  =  C 

...  ...  DO  40  K  .=..  1*J _ . _ _ _ _ _ _ _ _ 

. CHEC.K_EER-_.ZERC_RAIES _ _ _ _ _ 

_  I  F.LDWDT  (NW,K )..  LT  .T.  0  E-J.CX  GG-.TO  ,j6.5 . _ _ . . .  . . . . 

RATE  ( NW ,  K  )  =;  AL0G10 (  CW CT  ( NW  ,K  )  )  58 

_ RTEMP(NW.K)  =  1.0/(TPCLY(NW»K)  +273 .16) _ _ _ _ _ _ _ _ _ _ 

SUMXX  .=  .PARTIAL _SUM.  OF  X.SAUAR.E0  EIC. _ 

_ SUMXX  SUMXX ..+.  RTEMP  C  NW  ,K1M2 _ _ __ _ 

SUMYY  «  SUMYY  +  ( R ATE ( NW , K ) ) **2 

. .  SUMX  =  SLMX_3.  _RIEMP_(NW,)a _ _ _ _ 

SUMY  =  SUMY  +  RATE ( NW , K ) 

40  SUMXY  =  .SUMXY  RTEMP.I NW, K ) *R ATE( NW ,KJ _ 

GO  TC  55 

SET  UP  DUMMY  FCINTS  FCR  GP  I  FA  DWCT~~ VALUE  TwERO  ~ 

65  ACTE ( NW  )  =  0. 

_ ..P.REXI  NWX-.=  UU _ 

RATE  ( NU , K )  =  0. 

_ RXEMP.  (NV.K.)  _  =  .  a.OU15. . . . . _ . . . 

GO  TC  45 

_ 55  SLOP F  (NVl)  =  1 XJ4SUMXY- SUMX* SUMY  )/  (  X J*SUMXX-SUMX*<2  ) _ 

SPS(NW)  =( (SUMYY-(SUMY*SUMY/XJ)-( ( X J*$UMXY-SUMX*  SUMY )**2/ 

.  I X  J  *  XJ  *  SUM  X  X-XJ*SUMX*SU  KX.L  1  .ILL  XJ-2  .0)  1 _ __ _ _ _ _ _ 

ALPHA  =  ( SPS ( AW ) / ( SUMXX- ( SUMX*SUMX /XJ ) ) )*4 . 576 
- IE  (ALPHA)  58+58^51 _ _ 

57  SDS(NW)  =  SCRT ( ALPHA )  85 

- . -  _  GO  TC  55 _ _ _ _ _ _ - . . . . . . 

58  SDS(NW)  =  0.0 

59  8  E T A  .  =  l  S  PS (K W )* S  U MXX / ( X  J* S  U  M X  X-SU  MX  *SUMX  11. . . . 

IF(BETA)  62,62,61 

- 6i_SOT4-NW) — =— SORT  (SETA) _ : _ 94. 

GO  TC  63 

_ 62.  SO  I  (NW)  =  0.0 . __ . . . . . . . . _ . . . . . 

63  ACTE ( NW )  =  -SLCPE ( NW ) *4 . 576 

_ PREX  LN.W)  ._=„.  (  SLMXX4SUMY-SUMX1S.UMXY.)  /  (‘XJ* SUMXX- SUMX* *2) _ ! _ _ 

IF (NW. LT .20 )  C-C  TO  45 

- - IE-iNW.GT.60)  GC  TO  45 _ _ 
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TSUM  =  TSUM-SLGPE(NW) 

_ S£R£X  =.  S£R.5X  ♦  PBEMNWl . . . . . . . . ___ 

N  =  N+l 

_  45  CONTINUE _ _ _ _ _ _ _ _ _ _ _ 

C 

C  CALCULATE  AVERAGE  ACTIVATION  ENERGY  ANC  PRE-EXPCNENTIAL  FACTOR _ _ 

C 

AVPREX  «  SPREX  /  FLOAT ( N ) _ 

AVEA  =  TSUM/FLCAT ( N ) 

_ AVACTE  «  AVEA*  1.937 *2.303 _ _ _ _ _ _ _ 

C 

_C _ S_TART  LOOP  TO  CALCULATE  AFW _ _ _ _ _ 

C 

_ P_Q_LQ  _NW  .  _=_4j  9.8 _ _ _ _ _ _ . . .  .  . . . . . . . .  . . _____ 

Z  =  0 

_ DO  90  K__=  1»J  _ _  _ _ _ _ _ _ _ _ 

AFW(K)  *  RATE (NW*K )  +  AVEA*RTEMP(NW,K ) 

90  Z  =  Z  +  _AFW(K)_  _  __  _  _ _ _ 

FW(NW)  =r  Z/XJ 

_ WN.  ;  _f.LCAT.lNWJ _ _ ....  _ _ _ _ _  _ 

GG  =  ALCG10 ( 100 .-WN )  126 

SD  =  0  _  _  _ _ _ 

DO  93  K  *  1*J 

_93  SD  *  SD  ♦  (FW(NW)-AFW(K) )*»2 _ _ _ 

YK  =  J-l 

SPAFW  =  SORT (SO/YK) _ 

C 

C  WRITE  OUT  RESULTS 1  PERCENT  WT .  LOSS,  ACTIVATION  ENERGY,  PRE-EXPONENTIAL  _________ 

C  FACTOR,  AVERAGE  FW,  ANC  STANCARD  DEVIATIONS,  ALSC  LOG  WEIGHT  REMAI NINGI GG) 


X _ 133. 

70  WRITE  (6,1400)  NW, ACTE ( NW ) , SCS ( NW )  , PREX (N W ) , SD I ( NW ) , FW< NW ) , SDAF W , 

_ _ _.£LG _ _ _ _ _ 114 

WRITE  (6,1425)  AVACTE  141 

_ MJR.IXE. J  6,1415)  AV£R£.X _ _ 142 


WRITE  (6,1440) 

C  SETUP  INFORMATION FOR  CPSUEROUTINE,  SEE-OTHER  TrOGRA^S 

X _ _ _ _ _ _ litl 

L  =  3 

_ LS.  =  5 _ _ _ _ _ _ _ _ _ _ 

LW  *  100 

_ _ LN_*_  5_Q _ _ _ _ _ 

M  =  J 

DATA  A/  IF  .  /  _ _ _ 

JN  =  1 

.£___ . . . . . . . . . . . . . 

C  START  LCCP  FCR  PLOTTING  GRAPES  AT  10  PERCENT  WEIGHT  LOSS  INTERVALS 

DO  200  NW  =  1C, 99, 10 

_ _ _ DO  .100  X  =  1,  J _ _ _ _ _ _ _ _ _ _ _ _ _ 

X  ( K)  =  RTEMP(NW.K) 

_ ICQ  Jf  CXJ  . WjJIJ _ _ _ _ _ 

WRITE  (6,3000)  160 

_ _ _ _ WRIT  E_16.j_L.7QQ  !jNW _ _ _ _ _ _ _ _ _ _ _ 

C 

X. . . PLOT  GRAFH_C£.  LOG  .l.RAT O.F.  H E  IQ.HT  .LOSSl-AGA  IN.SI . RECJ.P.ROCAL _ _ _ _ — 
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C 

X. 

OF  TEMPERATURE 

161 

200  CALL  GP  (X,Y,L,LS,M,.JN,LW,LN,A,PLOT) 

M  =  100 

163 

WRITE  (6.3000) 

WRITE  (6.3100) 

167 

C 

C 

PLOT  GRAFH  OF  ACTIVATION  ENERGY  AGAINST  PERCENT  WT.  LOSS 

C 

CALL  GP  (AA.ACTE.L,LS.MiJN.LW.LN,A,PLOT) 

168 

169 

WRITE  (6,3000) 

WRITE  (6,3200) 

170 

C 

X 

C 

CALL  GP  (AA.PFEX,L,LS.M«JN«LW.LN,A,PLOT) 

171 

172 

WRITE  (6,3000) 

WRITE  (6,3300) 

173 

174 

00  75  1=1,87 

-  BB  ( I )  =  ALGGlQ  (-LQ0^-AA(  It3J  I 

179 

75  WF(I)  =  FW(I+3) 
_ LW  =  95 _ 


M  =  87 


X _ 


C 

PLOT  GRAFH  OF  LOG(AFW)  AGAINST  LOG ( PERCENT  RES I CLE  WEIGHT) 

X  - 

CALL  GP(EB,WF,L,LS,M,JN,LW,LN,A,PLCT) 

188 

WRITE  (6.3000) 

189 

LW  =  101 

M  =  101 

DATA  B/1H1 , 1H2 , 1H3, 1H4, 1H5, 1F6, 1H7, 1H8/ 

AA.1100)  =  0.0 _ 

AA(lOl)  *  100.0 

DO  -LQL— _ _ _ 

DWDT(IOO.K)  =  0.0 


1C1  DWQJllOl  ,  K )  =  ox 


WRITE  (6,3400) 

TALL  GP  (AA.DWCT.I.LS.M.J.LU.LN.R.  PLOT  ) 

199 

200 

X 

WRITE  (6,3000) 

C 

X  _ 

LOOK  FOR  FURTHER  SETS  OF  DATA 

201 

3C0 

READ  (5,1300)  MORE 

T  F  (  MORF  .  FQ  .  3  )  Kfl  Tn  1 

202 

80 

.  4X00 

STOP 

.  FORMAT  1 2X ,  A3 . 2  X . 8  A6 ) 

11C0 

L2C0 

FORMAT  (10X,A3,2X,8A6) 

FORMAT  (H  ,A5fl4fF13.«5,F6.1fF13.5,F6.1,F13.5tF6.1) 

1300 

1400 

FORMAT  (ID 

FORMAT  (1nX.I3,4X,-3PF7,3,5X,F6.3.5XTflPFA.3f7(5Y.F6.3)t2(5XrF6.411 

1425 

1435 

FORMAT  (//10X.29H  AVERAGE  ACTIVATION  ENERGY  =  .-3PF7.3) 

FORMAT  (10X.17H  AVFRAGF  LOG  PRFX.10X.JHa  .F6.3) 

1440 
_ 1500 

FORMAT (1CX.34HB0TH  FOR  20-60  PERCENT  WEIGHT  LOSS) 

XJQRMAT  (/8X,7HWT  1  OSS  , 2X.  8HFA  I  KO.AL  I  J3X  .‘8HST  .DF  VN  .  ,  3X  .  8H1  OG  PRPX, 

mm 

.3X.8HST.CEVN. .2X.10HAV.LOG  AFW.2X, 8HST .OEVN. , 2X , 1 1HL0G  RiES.WT.) 

-EQ&MAX  110X.lflHI.0G  RATE  VS  1/TEMP/ 1  OX.  14HWFTGHT  LOSS  *{.141 _ 
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1800  FORMAT  ( /  10X » 1 1HRUN  IC  NOS  »9(A5,2F.,l  )) 

1 9C0  FORMAT  ( 1  OX « 13HERR0R  FOR  W  «,I4.7HRUN  NO  .A3.6H  READ  ,I3< _ 

. 9H  I NSTEAD. ) 

2000  FO R  M AT_J_10 X j.2 5HLES S  THAN  3  HEATING  RATE S/1H1) _ 

3CC0  FORMAT  (1H1) 

3100  FORMAT  ( 10  X  >32  H  A  CTI VAT  I Q N  ENERGY  VS  WEIGHT  LOSS  > _ 

3200  FORMAT  (10X>22HPRE-EXP  VS  WEIGHT  LOSS) 

3300  FORMAT! IPX, 46FAVER  LOG  AF(W)  VS  LOG  PERCENT  WEIGHT  REMAINING) 

3400  FORMAT  110X »48HC0MPQUNDE0  RATE  OF  WT.  LOSS  VS.  PERCENT  WT.  LOSS  ) 

_ EM? _ _ _ 
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